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Abstract. Recent lattice studies of hadron properties, in particular of exotic states and charmonia are 
reviewed. Sea quark and quark mass effects are discussed as well as decays and mixing. 
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1 Introduction 

While the simplicity and elegance of QCD is very ap- 
pealing theoretically, the phenomenological observations 
of spontaneous chiral symmetry breaking and even more 
so of the confinement of colour charges turned it into 
a major calculational nightmare: it took almost twenty 
years after the discovery of asymptotic freedom to convinc- 
ingly demonstrate that the QCD Lagrangian indeed im- 
plies these highly non-trivial collective phenomena. This 
was done by numerical simulations; an analytic proof is 
still lacking. 

Fortunately, due to the property of asymptotic free- 
dom, many short distance/high energy QCD problems can 
be approached by means of perturbation theory. This need 
not be so since the very basis of the perturbative expan- 
sion is shaky: confinement implies that quark and gluon 
fields never appear as asymptotic states. Fortunately, the 
success of jet phenomenology suggests that QCD is rea- 
sonably benign in the high energy region. This is very 
different in the low energy regime of strong QCD. Among 
the few analytical tools that exist are the strong coupling 
and the expansions, effective field theories (EFTs) 
and various QCD inspired or phenomenological models. 

QCD can in principle be solved rigorously by means 
of lattice simulations on a computer. In practise how- 
ever computational resources are finite. This leads to pion 
masses of typically more than 400 MeV or to simulations 
within the quenched approximation, where sea quark ef- 
fects are neglected. While these systematic uncertainties 
reduce with faster computers, improved numerical algo- 
rithms and theoretical ingenuity, it will always remain de- 
sirable to combine lattice simulations with EFT methods 
or, where necessary, with QCD motivated models as only 
relatively simple questions can directly be addressed on 
the lattice. In many cases problems of phenomenological 
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interest factorise naturally into a high energy electro- weak 
part and a low energy QCD part which can then be evalu- 
ated on the lattice, for instance electromagnetic form fac- 
tors and weak decay matrix elements. In many cases QCD 
problems can also be factorised within the framework of 
EFTs into low and high energy parts. Lattice simulations 
turn out to be invaluable to gain insight into the dynamics 
of QCD as many parameters such as the number of active 
quark flavours, quark masses, number of colours, temper- 
ature etc. are not limited to their phenomenological values 
but can be varied. 

Several good reviews of lattice calculations exist in the 
literature and I refer to them for details on theoretical 
aspects, calculational methods and wider phenomenologi- 
cal implications . The topics covered in Chris Michael's 
review j^j "Exotics" have a non-vanishing overlap with 
this article, albeit written from a (slightly) different per- 
spective. In my talk I covered the baryon spectrum and 
structure as well. Due to the page limit (and lack of new 
content) I refer to my recent review Subjectively se- 
lected highlights since then include two studies of gener- 
alised parton distributions and a determination of 
the pion form factor j6j. Some progress has been made in 
the calculation of electromagnetic to Z\ transition form 
factors a new study of the mass spectrum of excited 
nucleons has appeared ^ and a first step towards consol- 
idation of previous results was performed |^. I will not 
mention the 0"*", because there are no lattice results. In- 
stead I will concentrate on the spectroscopy of charmonia 
and exotics. 



2 Glueballs and friends 

It was realised as early as 1974 ITU] that QCD offered the 
possibility of bound states composed only of energy. These 
closed-string-/glue-states/boxcitons became subsequently 
known as glueballs [111 . It is hard to imagine anything 
that demonstrates confinement more than the discovery 
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of quarkless massive bound states. One can also fanta- 
size that non-perturbative physics with particles similar 
to glueballs might play a role in future theories beyond 
the Standard Model. After all QCD is the only part of the 
Standard Model with a chance of a mathematically rigor- 
ous definition. So glueballs, the simplest possible colour- 
neutral states, are a challenge that certainly has to be 
addressed. 

The situation is complicated theoretically as well as 
experimentally by the possibility of mixing with standard 
quark model states. Even worse QCD does not know what 
is exotic and what is not. The very meaning of this term is 
"not quark model" and it is not at all a priori clear what 
exactly we mean by quark model. For instance we can ask 
ourselves how closely the quark model describes something 
as profane as the QCD proton. While introducing the con- 
cept of constituent quarks eliminates the puzzle that more 
than 98 % of its rest mass is generated by spontaneous 
chiral symmetry breaking, i.e. by the glue, it is impossi- 
ble to attribute the proton's spin and momentum exclu- 
sively to quark degrees of freedom. Gluodynamics plays 
an even greater role in properties of the rj' and the tt. So 
if the gluons already leave big footprints on quark model 
states how can we distinguish these from glue states? One 
possibility that immediately springs to mind would be to 
look for spin-exotic states: when coupling spin and angu- 
lar momentum within mesons, only certain combinations 
of J^*-^ are allowed while for instance 0^^, 0^ , 1 ^ , 2^ 
are forbidden. Again the situation is complicated by the 
possibility of qqqq molecular states to couple to these same 
quantum numbers. 

One toy model is the pure Yang Mills theory of gluo- 
dynamics. In this fictitious world no quark fields exist but 
many central features of QCD including confinement and 
asymptotic freedom are still reproduced. This well defined 
and self consistent theory yields a rich spectrum of glue- 
balls, all of which are absolutely stable, with the exception 
of very heavy ones that decay into lighter glueballs. 

One can go one step further towards QCD by includ- 
ing quark fields that propagate in the gluodynamic back- 
ground but whose feedback onto the vacuum is neglected, 
the so-called quenched approximation. Quenched QCD 
(qQCD) is no quantum field theory since unitarity is vio- 
lated, however, even in this approximation chiral symme- 
try is spontaneously broken. Neglecting the feedback of 
slowly moving heavy quark sources onto their environment 
is a very natural thing to do. Quenching is also justifiable 
in the limit of the number of colours N — > oo, however, it 
is not always clear whether 1/3 ^ 1. Simulations of pure 
SU{Nc) gauge theories for pitHI = 2, . . . , 6 seem to 
indicate this. As it should be ratios of light hadron masses 
from lattice simulations of qQCD have been found to be 
inconsistent with the observed spectrum however the 
differences are typically smaller than 10 % suggesting that 
the quenched approximation has some predictive power if 
cautiously consumed. The consequences of violating uni- 
tarity at light quark mass can become dramatic in some 
channels and in particular in the scalar sector |15| : roughly 
speaking as the axial anomaly does not exist in qQCD the 



rj' will be a surplus light Goldstone boson. The impact of 
this can be investigated in quenched chiral perturbation 
theory where diagrams that include transitions from scalar 
to two TTS yield unwanted contributions that explode as the 
TTS, including the would-be 77', become light. 

One justification for quenching is that the computa- 
tional effort is easily reduced by a factor of 10^. So it makes 
complete sense to learn and to understand statistics and 
systematics from a quenched study, prior to doing the real 
thing. Moreover, most models used in hadron physics ne- 
glect quark pair creation and annihilation, so the model 
builders can still learn from qQCD. Also, "un-quenching" 
a model to compare it with lattice results is easier than un- 
quenching the lattice simulation. Last but not least, life 
is easier in quenched as lattice studies of strong decays 
and excited states are notoriously complicated: glueballs 
will not mix with quark model mesons, spin exotic meson- 
gluon hybrids are distinct from mesons, are stable and do 
not mix with four-quark molecules either. 

2.1 Heavy glueballs and charmonia 

In Fig.^we compile the glueball spectrum (filled boxes) of 
gluodynamics [T^ltTC] . The scale Tq ^ !^ 394 is set from po- 
tential models. Since this is not the real world a systematic 
scale error of about 10 % should be assumed, for stable 
states that do not mix! While these newer results agree 
with the spectrum of Ref. (17) . during the past decade 
the statistical errors have been reduced by factors varying 
between 1.5 and 2, depending on the channel. It has also 
been established that the 3++ glueball is lighter than the 
1++. This had only been suspected in the earlier reference, 
where the J assignment of this state remained ambiguous. 

The lightest glueball is a scalar, followed by a tensor 
and a pseudoscalar. Possibly the lightest state with exotic 
quantum numbers is a 2"' glueball around 4 GeV. Quali- 
tative features can be understood in a bag model while the 
flux tube model seems inadequate ^S]. Without numeri- 
cal simulation it would have been completely impossible 
to come close even to a semi-quantitative understanding 
of these truly unconventional bound states. 

We include experimental charmonium levels (lines) for 
comparison. The notation refers to these states. 

We have only included confirmed resonances, with masses 
taken from the Particle Data Book. For the 77^ we took 
the Belle result from double charmonium production I19j , 
3630(8) MeV. One often thinks of the charm quark as 
heavy in the sense that its mass is much bigger than 
mesonic and baryonic QCD binding energies. However, 
the spectrum of states entirely made out of glue with no 
quarks at all covers a similar energy range! We also display 
recent quenched lattice results from two groups | i20u2m^ 
(CP-PACS and Columbia). These have been obtained us- 
ing relativistic charm quarks on anisotropic lattices with 
the ratio of spatial over temporal lattice spacings ^ 2 
(Columbia) and ^ « 3 (CP-PACS). A recent simulation 
of J < 1 5 and P wave charmonia on isotropic (C = 1) 
lattices by QCD-TARO j^Sl (who also study charmonium 
wave functions) confirms these findings. Note that in all 
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Fig. 1. The quenched charmonium spectrum (CP-PACS Columbia |21I22| '). glueballs |16I13| and spin-exotic cc-glue 

hybrids '22|, overlayed with the experimental spectrum. 



these simulations the effect of diagrams with disconnected 
quark lines has been neglected. One might expect OZl vi- 
olating contributions from these, in particular for states 
that lie close to glueballs with the same quantum numbers. 

The charm quark mass has been adjusted such that 
the spin averaged IS" state reproduces experiment. This 
means that in this case the quenching errors of up to 
10 % apply to spin-averaged level splittings with respect 
to this ground state, rather than with respect to zero en- 
ergy. This systematic scale uncertainty caused by omitting 
sea quarks renders it irrelevant whether we use a "con- 
structed" scale like tq or an experimental mass like the 
IP — IS* gap as an input (which would increase the split- 
tings by about 6 %). For the fine structure, potential mod- 
els combined with lattice results j21] tell us that we should 
expect to undershoot the real world number by up to 40 
%. For states above threshold or where mixing effects and 
strong decays play a big role radical changes might occur 
while for the 2P excitations finite volume effects could be 
an issue. 

There are different ways of obtaining a given JP^: for 
instance J^*^ = 1^^ can be an S* or a D wave. The corre- 
lation function associated with the respective D wave op- 



erator decays very fast in Euclidean time, into the same 
ground state as the respective S wave: it appears that 
the charm quark is too light to turn L into an (approxi- 
mately) good quantum number. A similar behaviour has 
been observed for the 2++ and states (either P oi F 
waves). 

In addition to the standard charmonia and glueball 
states the figure contains the lightest two spin exotic cc- 
gluon hybrids . At least two other studies of hybrid 
charmonia with relativistic charm quarks exist to-date |25l 
1^ . Given the fact that not even L is a good quantum 
number it is not clear how one would distinguish non spin- 
exotic hybrids from conventional radial excitations. The 
lightest exotic hybrid turns out to be a vector, 1 ^ . This 
is followed by 0"' and 2"' hybrids. 

Once sea quarks are switched on, most of the states 
calculated on the lattice will decay strongly and mixing 
will occur too. In this context it is interesting to see that 
the spin-exotic 0"' glueball and ccg hybrid have simi- 
lar masses. We have included two experimental thresholds 
into the figure: 1 charmonia can and will decay into a 
DD meson pair. Nonetheless resonances that exceed this 
threshold by almost 1 GeV have turned out to be experi- 
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mentally detectable. The exotic 1 ^ hybrid cannot decay 
into an identical boson-antiboson pair since PC = — . The 
next possibility, a D and a vector D* , is suppressed both in 
flux tube models |27j where the quark-antiquark pair that 
breaks the flux tube is produced in its centre, as well as 
in the heavy quark limit in which the angular momentum 
of the heavy quarks is fixed. The next possibility, which 
follows this S + P selection rule, would be a decay into a 
D and a P wave isovector D**[~ £>i(2420)]. Interestingly 
the mass of the 1 '"is compatible with this threshold such 
that mixing with (would be) molecular states is an issue. 
Should the state be below this threshold a resonance with 
a width smaller than 100 MeV might be possible, with 
dominant decay into a Xc under emission of a light scalar 
meson |2H1- 

At present we do not know whether the overlap with 
glueball states in the quenched approximation has any im- 
plications on the phenomenology of charmonia. Naively 
the phase space for decay of heavy giueballs appears huge 
with no quarks to preserve and more than 3 GeV of en- 
ergy to disperse. However, before drawing definite conclu- 
sions dynamical issues related to the glueball wave func- 
tions need to be investigated. In particular mixing is a 
strong possibility: for instance the 'tp{3770), il}{2S), DD 
threshold and the vector glueball all lie within the same 
100 MeV mass window, allowing for a potentially interest- 
ing phenomenology, in particular since there is little evi- 
dence from lattice simulations to support distinguishable 
25* and 11? levels in this region. Recently, the interesting 
question of a (non-exotic) hybrid component in the IS" 
states has been addressed in a very nice study "SS" , in the 
framework of NRQCD. The conclusion is that this contri- 
bution is weak. The question becomes more exciting but 
also tremendously less approachable if the 25" charmonium 
state and glueball channels were considered in addition. 

With data from Belle and Babar emerging, CLEO-c 
being online and the possibility of BES III starting to take 
data in 2006 the charmonium region is also exciting exper- 
imentally. Unfortunately, there is no way of producing say 
spin-exotic 1 states at any detectable rate in a decay 
starting from a vector resonance. For this we might have 
to wait for the proton-antiproton PANDA experiment at 
GSI. However, in the meantime there will definitely be 
progress in measuring spectrum and decay rates. Possibly 
the glue-richness of this mass region will leave its imprint, 
for instance by enhancing OZI suppressed processes. 

2.2 Light scalars: today 

A lot of experimental attention has been devoted to the 
spectroscopy of light scalar mesons. The reasons are three- 
fold: the scalar sector is intimately linked to chiral symme- 
try breaking which results in both light tts but at the same 
time also in heavy trs. Furthermore, the lightest glueball is 
predicted to be a scalar and, as an added extra, the light- 
est four-quark candidate state has scalar quantum num- 
bers as well: pseudoscalar-pseudoscalar bound states can 
be lighter then quark model scalars since chiral symme- 
try breaking makes the latter "artificially" heavy and the 



former "artificially" light. It is also not completely acci- 
dental that gluons play a dominant role in both the pseu- 
doscalar and the scalar sectors, in the former due to the 
axial anomaly and in the latter due to giueballs, both of 
which result in large OZI violations and splittings between 
/ = and / 7^ mesons. Unfortunately, it is exactly this 
rich phenomenology which complicates the interpretation 
of experimental results as well as theoretical calculations. 



The hghtest scalar "particle" is the a [or /o(600)]. It 
seems by now clear that this pole shifts the phase of the 
TTTT S wave, both from tttt scattering experiments as well 
as from proton-antiproton collisions with a Stt final state. 
With an imaginary part of the T matrix pole of 450(150) 
MeV, which almost exceedes its real part [800(400) MeV], 
few people are left who would call this sort of "resonance" 
a "particle" . However historically, before the advent of lat- 
tice simulations, the cr was among the first glueball can- 
didates. Next there is the /o(980) which in contrast has a 
width not much larger than 50 MeV, and is accompanied 
by the 1=1 ao(980) triplet. This "narrowness" is due 
to the fact that these states cannot decay into KK, be- 
ing 10-15 MeV lighter than this threshold. This proximity 
also turns them into natural candidates for either slls or 
KK bound states. For simplicity We will not distinguish 
between four-quark and meson-meson bound states but 
refer to both as "molecules" . 



Lattice results I17ll30lll6lll3j suggest that the lightest 
glueball should be a scalar with mass somewhere between 
1.4 and 1.8 GeV. While all raw lattice data agree within 
statistical errors of some 40 MeV, the large uncertainty 
quoted above is due to the scale uncertainty from using the 
quenched approximation. These results lie on top of three 
experimental scalar / = resonances, the broad /o(1370), 
the extremely well studied /o(1500) and the /o(1710). The 
standard picture is that these states are mixtures between 
a glueball and two / = nonet mesons, one with domi- 
nantly uu -\- dd and the other with ss quark content. The 
seven remaining / 7^ nonet members are most likely the 
four ^0(1430) as well as three ao(1450) states. Whether 
there is an excess of experimental states over quark model 
states or not critically depends on how resonances are or- 
ganised into nonets. It has to be said that the ao(1450) is 
quite broad and that the experimental evidence for these 
states is not rock hard. On the other hand, with a splitting 
of 450 MeV, it is extremely hard to reconcile the and 
the ao(980) into the same nonet, in particular the near 
degeneracy of the oq with the /o(980) suggests very tiny 
OZI violating effects. [This in turn also means that there 
is little chance for a significant glueball component in the 
/o(980).] The remaining puzzle is the absence oi S ^ Q 
partners of the ao//o(98G). Most likely these states have 
a large molecular component and are part of an inverted 
nonet [HJi together with the /o(980), a would-be- tttt a 
state and would-be- tt/C ks, which just happen to be ex- 
tremely fragile, due to the light tt involved. 
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2.3 Light scalars: future 

To support the above picture an unambiguous experimen- 
tal identification of the ao(1450) and a theoretical clarifi- 
cation of the molecular nature of the ao//o(980) states 
rank high on the wish list. Different experimental and 
theoretical inputs have led to various mixing models '32( 
13.3.1,34.1.35, • Amsler and Close suggest that the light 
quark state mainly goes into the /o(1370), with a sublead- 
ing component mixing into the /o(1710) while the ss and 
the glueball are mainly distributed between the /o(1500) 
and /o(1710). This view has been superceeded by Close 
and Kirk's more recent analysis of decays and produc- 
tion rates (HH] which suggests the /o(1500) to be domi- 
nantly glueball in character with light and strange quark 
components mixing destructively and additively into the 
/o(1370) and /o(1710), respectively. In contrast Refs. [HSl 
IM| predict the /o(1710) to be the (dominant) glueball- 
state with the two lighter resonances being composed pri- 
marily of the quark model mesons. 

All these predictions crucially depend on the input pa- 
rameters used, in particular on the ordering of "unmixed" 
states and on assumptions on the mixing matrix. One 
might expect the unmixed dd + uu state to have a mass 
close to that of the ao(1450) but unfortunately the lat- 
ter resonance is very broad experimentally and not overly 
well established. While in QCD with light sea quarks all 
states will automatically be "mixed" , the quenched ap- 
proximation is ideal to address the question of ordering 
of "unmixed" states and even for estimating off-diagonal 
mass-matrix elements. The main problem however is that 
in this case correlation functions will loose positivity at 
small quark masses as the tt and the rj' are degenerate P^ . 
In practise this means that at best ss scalars can be ap- 
proached using traditional methods. 

In Ref. significant finite size effects and a strong 
lattice spacing dependence are observed. After extrapola- 
tions the quenched ss scalar appears to be about 300 MeV 
lighter than the pure glueball. RBC PB' quote an even 
lower mass of 1.04(7) GeV for the light quark oq and 
0.9(1) GeV for the I = fo/a. The leading chiral correc- 
tion from the Trrj' loop was included into this study. One 
would however hope that systematic uncertainties will be 
identified and investigated in the future. Should this re- 
sult be confirmed, then the ao(980) could be explained 
as a quark model state but what is the /o(980) and why 
are the K^s, that are experimentally well established, so 
heavy? In a very sophisticated study Bardeen et al. [371 
I15j determine the relevant low energy parameters of the 
quenched chiral Lagrangian and remove the leading order 
quenched chiral power and log. After testing the approach 
in the pseudoscalar sector they are able to determine an 
ao mass of 1.33(5) GeV, expressed in units of the p mass, 
which is by about 200 ± 100 MeV lighter than the glue- 
ball, extrapolated to the continuum limit and converted 
into the same units. The unmixed ss state might then be 
expected to be fairly equal in mass to the glueball. The 
main difficulty is that with only two lattice spacings at 
hand the extrapolation to the continuum limit is still a 
bit shaky. 



The situation has also been studied by UKQCD [211 in 
QCD with sea quarks at lattice spacings a « 0.1 fm and 
a w 0.13 fm. In this case the oq is significantly lighter but 
the same holds true for the flavour singlet meson/glueball 
(see Fig. |21 below). 

In addition to the studies of quark model scalars some 
results exist on four quark molecules, but with degenerate 
masses, with the aim of addressing the possibility of a tttt 
bound state. Results of quenched studies PH] are incon- 
clusive due to the requirement of large volumes to avoid 
squeezing the tts on top of each other when they intend to 
be elsewhere: without a careful finite size study a residual 
interaction energy can hardly be disentangled from the tttt 
scattering phase |4()j . Diagrams with disconnected quark 
lines have also been neglected (with one exception j41|). 
In the absence of definite results in the quenched approx- 
imation, the SCALAR collaboration has started the very 
ambitious program of determining the mass of an / = tttt 
state incorporating sea quarks j42j . In view of the above 
it is highly attractive to study simplified cases, with heav- 
ier quarks and / = 1, like KK or DK (see also Sec. 12.41 
below) molecules, both with at least equally important 
phenomenological implications. 

Few exploratory studies of strong decays exist. Since 
we are working in Euclidean time there is no concept of 
asymptotic states. Neither can we calculate the imagi- 
nary part of a forward amplitude but there are bug-fixes: 
the method of choice is the computationally very chal- 
lenging exploitation of finite size effects introduced by 
Liischer (40) . The poor man's approximation has been de- 
veloped by Gottlieb et al. 43 and Michael jH^- In the 
latter case an on-shell transition matrix element is eval- 
uated by adjusting the rest mass of the decaying parti- 
cle to the energy of the outgoing state. In QCD with sea 
quarks this method can only be approximate as the opera- 
tor used to create the initial state will in general also have 
a non-vanishing overlap with the final state. The transi- 
tion matrix element is related to a coupling which in turn, 
assuming momentum independence, can be normalised to 
phase space, predicting the partial decay width. The on- 
shell condition implies that for a given mass of the final 
state the relative momentum has to be adjusted to guar- 
antee energy conservation. As momentum is discretized on 
a lattice this imposes some constraints. The initial state 
can be boosted such that at a fixed mass several relative 
momenta can be (approximately) realised, on sufficiently 
large volumes, and the momentum dependence of the cou- 
pling checked. Chiral perturbation theory can also help 
to connect results obtained at different quark mass, once 
these are sufficiently light. 

Promising results using this approach have been ob- 
tained recently for the p — > tttt decay in n/ — 2 QCD 
by McNeile and Michael jJH]. Decay couplings of a glue- 
ball to two pseudoscalars have so far only been computed 
by GFll PO] some 9 years ago: a mass dependence has 
been observed with a stronger coupling to heavier mesons, 
which then has to be folded with phase space. 

Elements of the glueball-scalar mass matrix have been 
calculated by two groups, again GFll in qQCD P3j and 
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Fig. 2. The scalar "glueball": qQCD vs. Uf ^ 2. 

UKQCD for rt/ = 2 ^H], for quark masses around the 
strange quark. Both collaborations obtain mixing energies 
of about 300 MeV and 500 MeV on relatively coarse lat- 
tices, respectively. Starting from two degenerate unmixed 
states this would imply level splittings between the mixed 
states of 0.5 - 1 GeV. An extrapolation to the continuum 
limit by GFll who also simulated three finer lattice spac- 
ings resulted in 61(45) MeV, a very reasonable value but 
with a 100 % uncertainty. 

We are still in the position that the combined "world 
data" on the scalar nj = 2 "glueball" , plotted as a func- 
tion of the squared lattice spacing a^, fits into Fig. 12 
The quenched case |13II17| is included for reference. The 
un-quenched results have been obtained by use of three 
different lattice discretizations of the Dirac action: stag- 
gered (HEMCGC 117] ), Wilson (SESAM gH) and clover 
(UKQCD gSIEHI)- The quarks are all not much lighter 
than the strange quark, the scalar meson is still stable 
and the wave function turns out to be very close to that of 
the quenched glueball 0HES1- Most Uf = 2 points clearly 
lie below the quenched line, however, there is certainly 
a slope in the results, such that the mass in the physi- 
cal a = limit does not contradict the quenched result. 
Within the SESAM data set there is an apparent discon- 
tinuity because different points have been obtained at dif- 
ferent quark masses; the "glueball" becomes lighter as the 
quark mass is reduced. Clearly additional studies at lighter 
quark masses and different lattice spacings are required. 



2.4 The LI+/2317) 

Recently a narrow Dfj{2317) state has been detected by 
BaBar 0^1) dominantly decaying into DsTt. This finding 
was confirmed by CLEO ^ and by Belle The lat- 
ter two collaborations also reported a narrow resonance 



around 2537 MeV, decaying into D* + tt. Both states lie 
by about 40 MeV below the respective DK and D*K 
thresholds. In contrast potential model calculations sug- 
gest [22] much heavier masses for the missing D^o and 
D'^^ P wave states, rendering these into broad resonances. 
These expectations have very recently received experi- 
mental support from Belle's observation of the missing 
D* and D[ states at 2308(17) (15) (28) MeV and at 2427 
(26) (20) (15) MeV, respectively |^. These states indeed 
strongly decay into Dn and Z?*7r, respectively, with widths 
of order 300 MeV. 

If the new Dgj mesons were P wave cs mesons, why 
are they so light but the corresponding cd mesons are not? 
Naturally this question invites speculation that the new 
scalar state might be of a DK (or 4-quark) molecular na- 
ture, somewhat resembling the /o/ao(980) system 135. Of 
course such ideas eventually have to be substantiated by 
a QCD calculation and indeed lattice results exist: sim- 
ulations in the static limit within the quenched approxi- 
mation and with sea quarks [SSI, simulations includ- 
ing NRQCD /HQET 1 /m corrections 'WT and simulations 
with a relativistic charm quark 58 ,59 . The interpretation 
of these results is controversial: I ^| observed that both 
effects, including sea quarks and including relativistic cor- 
rections to the static limit, increase the mass of a quark 
model cs P wave scalar state, pushing it above the DK 
threshold and into agreement with potential model pre- 
dictions. I then concluded that lattice results are incom- 
patible with a pure quark model nature of the new D^j 
states. On the other hand, while confirming an increase in 
the predicted mass when incorporating sea quarks, the au- 
thors of Ref. conclude that their results are consistent 
with the states observed by Babar and CLEO, however, 
within errors no disagreement is seen with the potential 
model predictions either. Interestingly, while 1/m correc- 
tions to the static limit are substantial for charm quarks 
and increase the 0+ - 0" splitting by about 25 % EH, the 
1+ — 1~ splitting is found to agree with the 0+ — 0~ split- 
ting within statistical and systematic errors of about 15 %, 
in quenched as well as with sea quarks [35], as suggested 
by chiral symmetry in the heavy quark limit j6()| . 

The same chiral symmetry argument would also apply 
to molecular states, such that the mere discovery of similar 
Bsj mesons would not help in discriminating between me- 
son and molecule. Dynamical issues need to be addressed 
and decays investigated. In particular electromagnetic de- 
cays, where the theoretical understanding is much better 
than for strong decays, would reveal information about the 
internal structure of these states. The narrowness of these 
resonances suggests that this need not be a completely 
hopeless enterprise. If the lightest Dgj state was domi- 
nantly a molecule, then a quark model J = resonance 
should exist in addition. The recent Belle [SSI evidence 
of the Dg(2308) suggests that another broad resonance, 
centred around 2.4 to 2.5 GeV with dominant decay into 
DK, might be detectable. 

Better lattice studies are needed to clarify these im- 
portant issues, and also with respect to the possibility of 
similar states in the B meson system. Of course at physical 
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Fig. 3. Splitting of hybrid meson masses with respect to the 
respective triplet S wave ground states for the light, strange, 
charm and bottom cases. X denotes n, K, D and B mesons. X* 
denotes the p,K* ,D* and B* vector mesons while XI stands 
for the (likely) 1+ states 6i(1235), A'i(1270), L>i(2400) and 
i3}(5732), respectively. 

light and strange quark masses lattice QCD should repro- 
duce the experimental spectrum. So to gain insight into 
the nature of these states a high precision quenched bench- 
mark study is required: in this approximation molecules 
and mesons are clearly distinct. In addition, the light sea 
quark mass dependence of both, a possible dls molecule 
and a cs meson has to be traced down into the region 
where mixing can set in. Fortunately, the DK system 
is more user friendly than its KK counterpart: heavier 
particles can be accommodated in smaller volumes and 
the binding energy that is suggested by phenomenology is 
comfortably large with about 40 MeV, rather than a mere 
10-15 MeV. It is also conceivable to calculate matrix ele- 
ments that are related to electromagnetic decay rates. 

2.5 More hybrids 

The spectrum of cc glue hybrids has already been dis- 
cussed in Sec. 12.11 above. Light hybrids have been stud- 
ied to some extent as well, with 61,62 and without [631 
EKll^ sea quarks. All these results yield the same or- 
dering as in the charmonium case (and in fact also the 
bottomonium case) with 1 being the lightest exotic, 
followed by 0"' and 2"' . The 1 is consistent with 
a mass of 1.9(2) MeV, with the corresponding strange 
quark exotic about 200 MeV heavier, in quenched as well 
as in un-quenched simulations. This is heavier than the 
7ri(1600) candidate. Mixing with molecular states is a pos- 
sible explanation and the feasibility of studying this has 
been demonstrated by MILC 26 . However, at present the 
quarks are still so heavy that for instance the combined 
mass of a bi and a tt is around 1.9 GeV too. Clearly the 



lattice calculations are incomplete. Molecules have to be 
included and the quark mass dependence of the mixing 
matrix has to be studied carefully. 

Heavy hybrids can be studied using NRQCD [HIII65II66I 
|67"68"29','18] or the Born-Oppenheimer approximation [691 
64,67,28, 12|. In the charmonium case simulations with 
relativistic quarks on anisotropic lattices have been pur- 
sued by several groups 22„23 , 25 as well as with isotropic 
lattices |26II62| . Even a simulation of bottomonium with 
relativistic quarks on an anisotropic lattice exists jZOl- 
With the statistical errors of present simulations differ- 
ences between quenched and un-quenched data cannot 
clearly be resolved. In FigureOlWe summarise the present 
estimates of the splitting of the 1 ^ hybrid with respect 
to the respective vector meson state: the flavour depen- 
dence is tiny. In addition three decay thresholds are dis- 
played. As detailed in Sec. 12.11 a strong decay into two 
pseudoscalars is forbidden while the decay into vector and 
pseudoscalar is suppressed, in particular in the heavy quark 
limit. Theoretically the hh hybrid is most clean-cut but ex- 
perimentally hard to produce. 

On the lattice electromagnetic matrix elements can be 
calculated |2;-{ll29j but strong decays are very challenging. 
To this end McNeile et al. JUj predict that if the light- 
est bottomonium hybrid is indeed below the B** B thresh- 
old the dominant decay channel should be deexcitation by 
emission of a scalar: H}, — s- x&tttt, with a width of about 
100 MeV. The same argument should also be valid in the 
charmonium case, where phase space would reduce the 
width even further, but not necessarily for light hybrids. 

3 Conclusions: from fiction to fact 

A combination of new theoretical methods and comput- 
ing technology has allowed us to arrive at the bound- 
ary between qualitative test of principle and quantita- 
tive prediction in the complicated area of flavour singlet 
physics, strong decays and mixing. A few years ago at least 
three major technical challenges had to be overcome: light 
quarks, sea quarks and disconnected quark lines. We are 
about half way through by now. With enough effort de- 
voted onto the topics covered in this article, quantitative 
predictions are possible within the time scale that is rel- 
evant for experiments like glueX at JLAB or PANDA at 
GSI. 
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